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Abstract 

The aim of this contribution is to review some aspects of theoret- 
ical studies of a family of quasi one-dimensional (QID) organic con- 
ductors known as the Bechgaard salts.ln order to make it personal, it 
will retrace the evolution of the author's interest in this field. 

The generalized aim was the calculation of the electrical conductiv- 
ity of the Bechgaard salts and gaining some knowledge on the equation 
of state and thermal properties of these materials. The review ends 
with some ideas about the possible future developement of the field. 
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1 Introduction 



I have entered the Institute of Physics (loP) in Beograd in the spring of 
1985. Several years before that I have graduated in physics at the University 
of Beograd, and after some time spent in looking (in vain) for a job in various 
research establishcmcnts was delighted by the fact that I was finally accepted 
in loP. I have entered the Laboratory for interdisciplinary studies, and my 
research field was broadly defined as " the study of the behaviour of materials 
under high pressure" . My immediate task was to calculate the electrical 
conductivity under high static pressure of several normal metals.If I suceed 
in the calculation,my job in loP would became stable. If not.... My first 
contact with the organic conductors dates from that hectic period, when I 
discovered the existence of these curious materials in the literature. 

My real encounter with the organic conductors occured in Paris in Prance 
a few years after that. In the spring of 1987., I was awarded a French governe- 
ment schoUarship for a visit of 8 months to Laboratoire des interactions 
moleculaires et hautes pressions (LIMHP) at Universite Paris-Nord.Thc aim 
of my visit was to learn experimental techniques of high pressure research. 

While working in LIMHP,I managed (through a friend of a friend from 
Belgrade) to arrange a visit to Laboratoire de Physique des Solides du CNRS 
on the campus of Universite Paris XI in Orsay.I was to be received by "some- 
one" named Denis Jerome. Even while I am writing this, after nearly 17 
years,I am smiling when I think that on the day when my visit was being 
fixed,this name did not mean anything special to me! However,the afternoon 
before I was due to go to Orsay,I received a reprint of a then recent pa- 
per from "Contemporary Physics" from which I learnt that the man whom 
I was going to meet actually discovered most of the facts related to quasi 
one-dimensional (QID) organic conductors. 

The rest of the story is " condensed history" . The following day I went to 
Orsay and met D.Jerome. He received me extremely kindly,we had a friendly 
conversation, focused on high pressure and organic conductors ,and some- 
where around noon we went to the canteen. While there, I started thinking 
that it would be "great" if I could somehow get the possibility to work for my 
Ph.D. in Jerome's group. However, I was in a dilemma how to ask him such a 
thing after meeting him for the first time in my life. Whether or not he felt my 
thoughts I do not know. The important thing happened when in the course 
of the lunch we came to delicious camembert cheese. 
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Jerome simply asked me whether I would be interested to come and work 
some time in his laboratory.! started laughing,explained him my toughts,and.. 
that was it. We agreed that I would come to Orsay and work for my Ph.D. there 
but that I would defend it in Belgrade. When my schollarship ended,! 
went back to Belgrade, completed and defended my M. Sc., prepared all the 
necessary administrative formalities, and in autumn of 1989.,! was in Paris 
again,with a schollarship for 12 months. 

Jerome's group was (and still is) one of the leading centers of stud- 
ies of the organic metals,in particular of the Bechgaard salts. The name 
of these materials is derived from the fact that they were synthetized by 
K. Bechgaard and his group in Denmark pp,and their generic chemical formula 
is {TMTSF)2X .The symbol {TMTSF)2 denotes a complicated molecule of 
the name of di-tetra-meta-tetra-selena-fulvalene, which is the basic ingredient 
of all the Bechgaard salts, and X denotes an anion such as P-Fg? FSO^, NO3.... 

When the Bechgaard salts were discovered, it was thought that they mar 
ked "the road" to be taken in the quest for high temperature superconduc- 
tivity. This turned out to be wrong, so the motivation for their study (both 
experimental and theoretical) somewhat changed. The main motivation for 
their study nowadays is the idea that they are the simplest possible cases 
(and therefore easier to understand) of correlated electron systems, which (it 
is thought worldwide) the high materials are. 

Time went quickly and ! worked on the experiment aimed at measur- 
ing the electrical conductivity of the Bechgaard salt {T MT S F)2F SO^ j2]. 
When the experiment was almost complete, and we were starting to plan 
the continuation of it with the introduction of impurities, a critical part of 
the experimental setup exploded. No one was harmed, but waiting for the re- 
placement ! started trying to do some theoretical work on the conductivity 
of the Bechgaard salts. Jerome noted it and put me into contact with !!einz 
J.Schulz of the theory group. A couple of months before that,on advice of 
Jerome ! have asked for the prolongation of my schollarship for another year. 
This was at first rejected,but the Laboratory (and Jerome personally) used 
their " contacts" in the administration. After the summer holidays,when ! got 
back to Paris,! found a hand written note from Jerome saying "Dear Vladan, 
Wellcome for another year in Paris".! was proud and happy,passed to the 
theory group to work with Schulz,and that is where my theoretical work on 
the organic conductors begun. 
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2 The calculation of the conductivity 



2.1 The first attempt 

Heinz Schulz was a tall, quiet, friendly man spending all the day in the 
laboratory,except when he had courses at the Ecole Politechnique.He was 
from Hamburg,got his Ph.D. there and came to Orsay as a post-doc. He pre- 
pared another Ph.D. there, and was immediately offered a job.Sadly,he died 
of cancer in 1996., and only then I discovered that we were of nearly the same 
age. We agreed that after my experimental work in Jerome's group it was 
logical that I continue by trying to calculate the electrical conductivity of 
the Bechgaard salts. 

Two immediate questions to which we had to find answers were the choice 
of the method of calculation, and the choice of the Hamiltonian (i.e., a theo- 
retical model of the material) which could be an appropriate description of 
the Bechgaard salts. 

The logical choice of the theoretical model of these materials was the 
famous Hubbard model. It was proposed several decades before the start of 
my work in Orsay,and it contains the basic "ingredients" for the correct 
theoretical description of a solid: atoms, localized in the nodes of some form of 
a crystal lattice and electrons localized on them jS] .Although it is intuitively 
simple and acceptable, to this day the Hubbard model was solved only in the 
simplest (one-dimensional) case P^.In one spatial dimension the Hamiltonian 
of the Hubbard model has the following form: 

N 

H =-tY. i.4+l,aCi,a + cl^Ci+i^a) + U ^h^'^Ui = Hq + Hi (1) 
i,(T=0 /=0 

The first term is the kinetic energy term (ifo))Cind the second one is the 
interaction term {Hi). All the symbols in this equation have their standard 
meanings: t is the electron hopping energy,the terms in parenthesis in Hq 
are the creation and anihillation operators for electrons on a lattice site i 
with spin a.In the interaction term,?/ is the interaction energy of a pair of 
electrons having oposite spins within an atom on lattice site I. 

A starting idea for the calculation of the conductivity was to try and apply 
the historic Kubo formula 0.1n modern transcription, such as Kubo's 
formalism gives the following expression for the real part of the conductivity 
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tensor: 

and / denotes the following 

/oo 
exp^'- {M^{t)MM)o (3) 
-oo 

where is the time dependent dipole moment operator.In the Heisenberg 
picture,this time dependence is expressed as 

M^{t) = ex/^°/^M^exp-'*^°/^ (4) 
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To complicate things further,the calculation of the correlation function in 
the expression for I demands the knowledge of the partition function and the 
density matrix of the system. Obviously,the Kubo formula although being 
exact is too abstract for performing a realistic calculation. 

2.2 The memory function approach 

After we jointly concluded that trying to apply Kubo's formula was actually 
leding to an impasse Schulz suggested that I try applying the so called "mem- 
ory function" method. Of course I accepted, and had an immediate feeling that 
I was entering a sort of scientific "terra incognita". 

The memory function method has its roots in work in pure rigorous statis- 
tical mechanics around the middle of the last century (|7|,[H| and earlier work 
cited there). A theory of many-particle systems was developed in these papers 
with the aim of describing transport properties, collective motion and Brown- 
ian motion from a statistical-mechanical point of view. This formulation, taken 
from the abstract of j7j, does not sound complicated. However, retracing in de- 
tail the calculation performed by Mori is a task of considerable mathematical 
complexity. The main point of using the memory function in any calculation 
of transport properties is that it gives the possibihty of evading the problem 
of formulating and solving the transport equations. Instead of dealing with 
the transport equations, this method gives the possibility of expressing the 
conductivity in terms of a regular " memory function" . 

Papers [7] and jH] are, in some sense, comparable to j3]:they are exact in 
principle, but apphcable to real calculations with extreme difficulty. However, 
some time after the publications of Mori, a technique was worked out which 
gave the possibility of applying the memory function method to real calcu- 
lations ii,iini- 

The equations needed for the calculation of the conductivity within the 
memory function method are [TU] 



Xab{,z) =« A- B »= -i / exp'-"^^ < [A{t), B{0)] > dt 



(5) 



and 
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In these equations Up = ATm^e^ /m^ is the square of the plasma frequency, 
rZe, e and me are the electron number density, charge and mass,while de- 
notes the ratio ng/me, which is the zero frequency limit of the dynamical 
susceptibility xi!^)- 

The integral in eq.(5) is a general definition of the linear response of an 
operator A to a perturbing operator B. It is an analytic function for all 
non-real frequencies z llOj.The symbol A{t) denotes the Heisenberg repre- 
sentation of the operator A.Eq.(5) can "evolve" from a general definition of 
a response function into the definition of a current- current correlation func- 
tion by introducing A = B = [j, if]. The current operator is denoted by j and 
H is the Hamiltonian of the system under consideration. 

It was decided, in line with |Tr|, to use the Hamiltonian of the Hubbard 
model in the calculation. The Hamiltonian of this model |S] in one dimension 
is given in eq.(l). The current operator has the form 

3 = I](ctaQ+l,^ - ci+i^aCl,a) (7) 

The final step which needs to be made before the actual calculation of 
the conductivity of the Bechgaard salts is the determination of the chemical 
potential of the electron gas on a ID lattice.lt may not be obvious at first 
why is the chemical potential so important for this work. The usual research 
practice is that pure Bechgaard salts are described as systems with = O.Any 
possible doping (and experiments with doped specimen are very interesting) 
is theoretically described by deviating ^ from O.The chemical potential is 
a function of the band-filling ?T,,the inverse temperature (3 and the electron 
hopping energy t.A determination of an analytical expression for this function 
demands a choice of a theoretical model for the description of the electrons 
in the Bechgaard salts: Fermi or Luttinger liquid? 

2.2.1 Fermi or Luttinger liquid? 

The basic ideas of what is now (for several decades already) called the Fermi- 
liquid theory are present in physics for a long time - it can be safely stated 
since the time of Sommerfeld. He showed that various experimental data 
on metals (for example the low temperature behaviour of their specific heat 
or the electrical conductivity) can be understood by assuming that the elec- 
trons in a metal behave like a gas of non- interacting Fermions. Shortly after 
Sommerfeld, came the results of Pauli,Bloch and Wigner. 
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They showed that the paramagnetic susceptibihty of the non-interacting 
electrons is temperature independent, and that the interaction energies of the 
electrons at metallic densities are comparable to their kinetic energy.The 
"final touch" was given by Landau ^2] -In these works Landau intro- 
duced a new way of thinking about interacting systems, which became crucial 
for the developement of condensed matter physics. For example,the notions 
of quasiparticles and elementary excitations,were introduced in these pa- 
pers. Methodologically,they were extremely important, because in them Lan- 
dau introduced the idea of "asking useful questions about the low energy 
excitations of the system based on concepts of symmetry, without worrying 
about the myriad of unnecessary details" jl4j . 

When does the Fermi liquid theory break down? Physically,the answer is 
logical: it happens whenever one of the measurable quantities which it aims 
to calculate diverges ^3] .It has been shown in field theory that the Fermi 
liquid concept breaks down in one dimension. Technically speaking, some 
vertices which are in the Fermi liquid theory assumed finite, diverge in one 
dimension because of the Peierls effect (for example jTHD-It is also known 
that the excitations in ID are not quasi particles,but collective charge and 
spin degrees of freedom,each of which propagates with different velocity • 

Optical and photoemission experiments on the Bechgaard salts show de- 
viations from the Fermi liquid behaviour jTTj. However, the energy scales 
on which the deviations occur in these two groups of experimental data are 
different. Interesting conclusions have been reached from c-axis conductivity 
measurements of the Bechgaard salts ISj.It was shown there that the ID Lut- 
tinger liquid description of the Bechgaard salts breaks down below a pressure 
dependent value of the temperature, and that the Fermi liquid description is 
restored at low temperatures of the order of 10 K. 

This subsection is certainly not a complete review of the problem " Fermi 
vs.Luttinger".It's aim is to try to justify to the reader the choice which was 
made in this work: to use the Fermi liquid theory. The logical question is 
"Why?". 

There are several motives for this decision. First and simplest - chronologi- 
cally,this calculation of the conductivity of the Bechgaard salts was completed 
in 1997. and started while the present author still was in Orsay. Second, at the 
time when this calculation was being performed many experimental data 
which are now available simply did not exist. 
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This implied that it is wiser to perform the calculation within the Fermi 
liquid theory as a very well known theoretical framework. Experimentally,it 
was known that the Bechgaard salts are not strictly ID but quasi ID - the 
conductivity along the so called c-axis is not exactly equal to zero. There is 
also the aposteriori justification: results for the temperature dependence of 
the resistivity are qualitatively similar to experimental data. 



3 The calculation-the practical part 
3.1 The chemical potential 

It is known since |1] that the chemical potential of the electron gas is zero 
for the single-band half- filled Hubbard model. The band filling is defined as 
the rato of the number of electrons present per lattice site to the maximal 
possible number of electrons per site A^. This value is obviously equal to 2 
since two electrons with different spins can be present on the same lattice 
site. 

" = ^ (8) 

Theoretically speaking,the value of rie is a parameter which can be varied 
at will. From the experimental viewpoint, this value can be changed by doping 
the specimen with electron donors or acceptors. As the chemical potential 
enters the expression for the Fermi fermi function,the first practical step 
in calculating th electrical conductivity had to be the determination of the 
chemical potential of the electron gas on a ID lattice 

The starting point of this calculation is the following equation: 



1 

n 



nh J 1 + exp/3(efc — /i) 

where all the symbols have their usual meaning,p = hk and n is the number of 
electrons per site. The integration is performed within the limits ±iTh/s and 
s is the lattice constant. In the ID Hubbard model tk = —2t cos(A;s), which 
leads to 

dp=Mk=-{Ae-e^Y^'^de (10) 
s 
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Inserting eq.(lO) into eq.(9) gives 



n^ — r 1 -(4^2 - e^Y^'He (11) 

TTS J-2t 1 + expP{ek — jJi) 

The initial problem of determining the function = ij,{(3,t,n) has thus 
been reduced to the problem of solving the integral equation eq.( 11) .Instead 
of proceeding with the full rigour of the theory of integral equations,one can 
tackle the problem by using a suitable power series representation of the 
Fermi function. 

The Fermi function can be represented as 



-1 CO 

l + - - ^) - - .) (12) 

where the symbol A2k+i denotes 

_ 2(-l)^+V^+^(2^^+^ - 1) ^ 
A2k+i (WT2)\ ^^^^ 

is the step function, 5^*^+^ (e — fi) are the derivatives of the 5 function and 
-82/0+2 are the Bernoulli numbers. 
Inserting eq.(12) into eq.(ll) gives 



n 



-I 

ns J- 



e(/. - 6) - E A2k+iP-'^'^'S''^'{e - /x) 

k=0 



(4^2-62) -1/2^, (14) 



Eq.(14) is integrablc by using the following relation, valid for any function 
f{x) and its n-th order derivative: 



(15) 



Inserting eq.(15) in eq.(14) and limiting the summation to terms with 
k <2 leads to an equation of second degree in // which can be solved to give 
the following result for the chemical potential of the electron gas on a ID 
lattice: 



{I3tf{ns - 1) \t\ 



1.1029 + .1694(/3i)2 + M54{(3ty 



(16) 
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Obviously,for /im„ ,,^i/i=0, which means that this result has as its hmiting 
value the well known result of Lieb and Wu 

The reader may at this point be tempted to ask for an explanation of the 
limitation to terms having A; < 2 in the calculation leading to eq.(16).The 
explanation is that the idea being to obtain an analytical expression for the 
function t, n) going to terms with k > 2 leads to expression which are 
too complicated for any applications in further work. 



3.2 The electrical conductivity 

With the chemical potential being known, it finally became possible to tackle 
the problem of the electrical conductivity of the Bechgaard salts [201 • The 
Hubbard hamiltonian and the current operator are given by eqs.(l) and 
(7). The functions x{^) cr(-2) are both expressible in complex form as 
a{z) = ctr^z) +iai{z) and xi^) = Xr{^) oiis further assumes that 

the frequency is a complex function,with z = zi + iz2 and that Z2 = olz\ with 
a >- 0,it follows from eq.(6) that 



AziTfil + ia) 



Xr{z) + iXi{z) 
Xo 



(17) 



This expression can be separated into the real and immaginary part. Taking 
the special case a = (justificed by the fact that the frequency measured in 
experiments is real) gives: 

and 



aj = ^il-^\ (19) 

The ensuing steps in the calculation of the electrical conductivity are in 
principle straightforward. The expression "in principle" is here amply justi- 
fied,because inserting all the '"sub-results'" into the expression for xab leads 
to an almost intractable expression. 

The "practical" calculation (practical is here employed in the sense lead- 
ing to the function an) starts by determining the current-current correlation 
function. This can be obtained by inserting A = B = [j, H] into eq.(5). 
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Evaluating this commutator is facilitated to some extent by the decom- 
position of the the Hamiltonian of the Hubbard model indicated in eq.(l). 
Relatively easily it can be shown that [j, Hq] = and that 

(20) 

All the symbols have their standard meanings. Transition to k space can be 
performed by relations of the following form [21] 

cl = N-'/'Y.^xp{zhls)ct,. (21) 

kl 

In this expression is the number of lattice sites, s is th elattice constant 
and L = Ns is the length of the specimen. The temporal evolution of vari- 
ous operators can be introduced in the calculations by relation of the form 
Ck{t) = exp —ie{k)tck where (for the one dimensional Hubbard model) e{k) = 
—2)f:cos(A;s). Inserting all this into the expression for xab gives the following 
expression for the susceptibility xi.^) 

X{z)= mi{l/[{l + eMP{-^-2tcos{gm 

P,9,k,q 

(1 + exp(/5(-/i - 2tcos{k))))] - 1/(1 + exp{p{-fi - 2tcos{p)))) 
1/(1 + exp(/5(-;U - 2tcos{q))))]{Uty{azi + i{zi + 2t{cos{q) + cos{p) 
—cos{g) — cos{k)))){cos{p + g)/2){cos{{q + k)/2)[cosh{g — p) — 1]/ 

{N^iiazif + {zi + 2t{cos{q) + cos{p) - cos{g) - cos{k)))Y)) (22) 

In order to simplify the calculation at least to some extent, the summation 
was limited to the first Brillouin zone, and the lattice constant s was set as 
s = l.This obviously implies that L = A^.The sum in the equation above were 
calculated under the condition a ^ because this condition is built in in the 
definition of x(-2^)-0nce the sums were calculated,taking into account that in 
reality the frequency is a real quantity,the limit a ^ was imposed. After 
performing all the summations and taking a — > the following approximation 
for the dynamical susceptibility was obtained. 
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X = (32i(-l + cosh{l))[l + exp{/3{-fi + 2t cos(l - Tr))))-^ 
(1 + exp{(3{-i2 - 2t))Y\UtlN^f cos2((l - 27r)/2) 
{Z2 + i{zi + 2t(-2 - 2cos(l - v:))))/{zl + {zi + 2t(-2 - 2cos(l - 7r)))=^ 

+ « 2267 » (23) 

The number << 2267 >> in the preceeding equation denotes the number 
of omitted terms. Obviously,such a long equation is inapphcable and has to 
be truncated after a certain number of terms. Taking the first 32 terms of 
eq. (23), multiplying out the products and powers, expressing the result as a 
sum,gives the real part of the dynamical susceptibility as 

Xr{z) = [128UH^cos''{{l - 27r)/2)]/[l + exp(/3(-^ - 2t)))''N\zl + 

+ {zi + 2btf)]+ « 527 » (24) 

and b = -4(1 + cos(l - tt)). 

Taking the first 20 terms of eq. (24) and imposing the condition lim Z2 ^ 
it follows that the real component of the function xr has the form 

XhM-Ej^ + I:,^ (25) 

The functions Ki and Lj can be read-off from eq. (24) developed to a given 
number of terms. Using this result, the immaginary part of the dynamical 
susceptibility (xi) is given by 

Xiiuj) = -2—P / — 26 

TT J [UJ^ — Uo) 

Inserting eq.(25) into eq.(26) and imposing the constraint i,j < 4 in eq.(25) 
leads to the following expression for xi 

Xi = {2ht/Tx){Ut/N'^f[ujQ/{ujQ + 2ht){ujl - (26^)^)] [4.53316(1 + 
exp(^(-;u - 2t)))-2 + 24.6448(1 + exp(/3(-/i + 2tcos{l - Tr)))))-^] + 
(2/7r)[cuo/(cu^ - {2htf)]{Ut/N^f x [42.49916(1 + exp(/5(-/i - 2^)))-^ + 

78.2557(1 + exp(/?(-/i + 2t cos(l - Tr))))-^] 

(27) 
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The final expression for the electrical conductivity follows by inserting 
this result into eq.( 18). After some algebra,one gets the final result for the 
electrical conductivity of QID organic metals 

= {l/2xo){u;lM[u;l - {htf]-\Ut/N^f x S (28) 

and the symbol S denotes the following function 

S = 42.49916 X (1 + exp(/5(-/i - 2^)))"^ + 78.2557 x 
(1 + exp(/3(-/x + 2t cos(l + 7r))))-2 + (ht/{ojQ + ht)) x 
(4.53316 X (1 + exp(/5(-/x - 2f))Y^ + 
24.6448(1 + exp(/?(-/x + 2t cos(l + Tr)))))-^) 

(29) 

3.3 Discussion 

Equations (28) and (29) are the final result for th electrical conductiv- 
ity of the QID organic metals obtained within the memory function the- 
ory.From the purely mathematical point of view they are "moderately com- 
phcated" ,the main " 'complicating factor" being the length of the two expres- 
sions. However, for applications in physics,the relevant question is whether or 
not the results of these two expressions agree with experiments. 

The main externally controllable parameters in experiments on organic 
conductors are the temperature, the doping and the frequency,and in the re- 
mainder of this chapter, the comparisons of the results of eqs.(28) and (29) 
will be performed with respect to these parameters. Model parameters which 
are constant in these two expressions were chosen as foUowsiA'^ = 150, C/ = 
At,ujp = 3U,Xo = 1/3 and ujq > O.GU.The lower limit for ujq was determined 
by imposing the condition that aji > 0,while all the other values were chosen 
in analogy with high temperature superconductors. The influence of doping 
(and its changes) on the conductivity can be studied through variations of 
the chemical potential which depends (among other parameters) on the band 
filling. A half-filled band (n = 1) describes a chemically pure specimen. Posi- 
tive deviations of the filling from the value n = 1 experimentally correspond 
to doping with electron donors. Negative deviations describe the doping of a 
specimen by electron acceptors. 
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The first test of tlie result obtained for the electrical conductivity was its 
application to the case of a ID half-filled Hubbard model. Namely,it is known 
from rigorous theory that its conductivity is zero (for example |2 25). Inserting 
n = 1 into eq.(29) and developing it in t as a small parameter, it follows that 
the conductivity is approximately given by 

aniuj) ^ 10-^(— )V(4.56 + .3(3t) (30) 

UJ 

Clearly,the numerical value of this expression is close to zero for physically 
acceptable values of the input parameters,of course excluding the case u = 
0. Looking from the experimental side,this can be interpreted as implying 
that weakly conducting phases of the Bechgaard salts can be described by a 
ID Hubbard model with a small deviation of the band filling from the value 
of 1/2. 

The experimental parameter which can be most easily controlled in exper- 
iments on organic conductors is the temperature of the specimen. The general 
conclusion of all such experiments is that the temperature dependence of the 
conductvity of the Bechgaard salts is extremely complex, and that it can not 
be explained by conventional theory of conductivity of metals. 

In the calculations discussed in this chapter no attempt was made to 
reproduce the experimental conductivity of any particular salt. However, the 
idea was to determine whether or not the one band ID Hubbard model can 
reproduce semiquantitatively the experimental data. Examples of real exper- 
imental data and theoretical calculations of the conductivity are available in 
[T7j,|2ni ^^cl many other publications. 

Equations (28) and (29) can be re-written as follows 

aniujo) = A [u^ - (6t)'] [Q + btZ{uJo + bt)-^] (31) 

where the sumbols A,Q and Z denote various frequency independent functions 
which occur in the expression for the conductivity.Developing this result in 
ujq up to second order terms,it follows that 



aRiuJo) = -A{Q + Z){bt)-' + AZ{bt)-'{uJo)-A{2Z-Q){bt)-\uJoy + ... (32) 

Fitting this equation to measured frequency profiles of the conductivity,one 
could determine the functions A,Z,Q and t. An approximate value of the 
static limit of the function aji{uJo) follows from eq.(32) as 
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an{u;o = 0)^-A{Q + Z){bt)~' 



(33) 



Fitting these two expressions to experimental data, it would be possible 
to determine the values of the functions A,Z,Q and the hopping integral 
t. Fixing all the parameters and varying the band filling would give the possi- 
bility of investigating the effects of doping on the conductivity.A preliminary 
investigation of this sort has been performed in |2Uj . 



3.4 The equation of state and specific heat 

It is hoped that the reader has by now gained an idea about the complexity 
of theoretical studies of the Bechgaard salts. The aim of this sub-section is 
to outline a determination of the specific heat per particle under constant 
volume of a degenerate electron gas on a ID lattice. Such a calculation may 
seem as an exercise in statistical physics and mathematics.lt fact, in recent 
years measurements of specific heat and thermal transport properties of high 
Tc superconductors and organic conductors have become a useful tool in 
studies of these systems. The hterature in the field is steadily growing,but 
; [21] ; are some useful examples. 

It will be assumed in the followin that the number of particles in the 
system is not conserved. Mathematically speaking,the starting point of the 
calculations is the equation 

dG = -SdT + VdP + /ic/iV (34) 

where G = fiN and all the symbols have their standard meanings. Inserting 
the definition of the thermodynamic potential G into eq.(34) and differenti- 
ating with respect to T,one obtains 

In the last expression the obvious transformation dP = ^dT (and similar 
relations for other variables) was applied. Eq.(35) is expressed in terms of 
the bulk parameters of the system. In order to re- write it in terms of the local 
variables, a change of variables —>■ nV and S —>■ nsV is necessary. After 
some algebra,it can be shown that 
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nV —(s ^ ii) ^ n^— 36 

We have here derived the differential equation of state (EOS) of any 
material.lt is at first sight complicated,but in applications to the Bechgaard 
salts it can be considerably simplified. Experiments on these materials are 
almost always performed under fixed volume conditions [2n]-This implies that 
all terms in eq.(36) containing volume derivatives can be disregarded. The 
final form of the EOS of QID organic metals thus emerges as 

dP , dn. , da. ds 

If = <« + sr'(-' + af' + "af 

The specific heat per particle is given by [2Zj 

-V = -{gf,)v - r(^). (38) 

Applying this definition and performing all the necessary algebra,one fi- 
naly gets the expression for the specific heat per particle 

T d^..dn d'^n T dn ds d'^fi 

n^^^) (39) 

In order to apply eq.(39) to the degenerate electron gas on a ID lattice, it 
is necessary to introduce into it appropriate expressions for the chemical 
potential, number density and entropy per particle. The chemical potential is 
given by eq.(16).As a first approximation,known results for the entropy per 
particle and the number density can be used [2H1-The entropy per particle is 
given by 

_ Q dF,,MT) 
n dT 
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The symbol F3/2{fi/T) denotes a special case of a Fermi-Dirac integral 
PHj.n is the electron number density and Q is a combination of known con- 
stants (such as the electron mass and the Planck constant). Finally,the num- 
ber density of a degenerate electron gas at low temperature can be expressed 
as 



n 



A X [l + S X T3/2 + ...1 (41) 



where the symbols A and B denote combinations of known constants. Inser- 
ting eqs.(16),(40) and (41) into eq.(39) gives the final result for the specific 
heat of the electron gas on a ID lattice. After all the necessary algebra, it 
would turn out to be highly non-linear, in line with the results of various 
experiments (such as Pl], (221 and later work). 



4 Conclusion 

The aim of this chapter was to review some results of theoretical studies of 
a class of organic conductors called the Bechgaard salts. From the start it 
was prepared as a "personal review" , which simply means that it is to a large 
extent based on previous research results of the author. At the end of such a 
review,one can reflect for a moment on the possible future developement of 
research on the Bechgaard salts, and the organic conductors in general. Like 
Lord Kelvin in the XIX century,one may be tempted to think that after 
more than 20 years since their discovery, most things about the Bechgaard 
salts are well known and only "small clouds" remain to be clarified.lt seems 
(at least to the present author) that this would be totally wrong,and many 
interesting problems are waiting to be solved. A few of them are mentioned 
in the follow up. 

The general question can be formulated as "Why should anybody on this 
planet study the organic conductors (and the Bechgaard salts as a paritular 
family of the organic conductors)?" The answer is "simple" : the Bechgaard 
salts are a QID example of strongly correlated systems, and therefore should 
be easier to understand than their 3D analogs, such as the high supercon- 
ductors. However, in spite of all the efforts,the physical mechanism of high 
superconductivity has not yet been discovered. 

A basic problem is,of course,the apphcability (or inapplicability) of the 
Fermi liquid model to these materials. Indications from various experiments 
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(for example [T7j) seem to be that the Fermi hquid is not a good descrip- 
tion of the organic conductors. To make things more interesting,there are 
also indications (such as j^Dlor JHl) that reasonable agreement between cal- 
culations performed within the Fermi liquid model and experiments can be 
achieved in some cases.A definite solution of this dilemma is certainly one of 
the problems waiting to be solved as soon as possible. 

The calculations reported in the present chapter were performed within 
the "original version" of the memory function method. The method itself has 
been considerably reformulated and modernized in the mean time, bringing 
it more "in line" with contemporary field theory. For a recent application of 
a modernized version of the memory function method see, for example. |3Uj. 
Another method which may be useful for calculations like those described in 
this chapter(but also for the determination of the phase diagram) has been 
developed in [31j. 

The determination of the phase diagram of the Bechgaard salts is an 
interesting problem on its own. The temperature dependence of their con- 
ductivity can vary in some regions of the P — T plane; accordingly,they can 
be insulators, superconductors but also normal conductors or semiconduc- 
tors. The general form of the phase diagram of the Bechgaard salts is known 
for several decades. However, much more complicated is the possibility of de- 
termining the parameters (such as t) of a Bechgaard salt from the analysis 
of the phase boundaries on a P — T diagram. For a recent example of such a 
work see [S2]- 

Interesting considerations have recently been made concerning the dimen- 
sionality of the Bechgaard salts. Namely,they are usually considered as one 
dimensional,or quasi one-dimensional. However, recent angular magnetoresis- 
tance osicllation (AMR) experiments [HS] on {TMTSF)2FS0^ have shown 
that a cylindrical Fermi surface can be formed for this material. Does this per- 
haps mean that the Bechgaard salts should be considered as Q2D materials 
is (to the knowledge of the present author) a completely open question. 

This list of selected problems concerning the Bechgaard salts could be 
continued,but it is hoped that the examples present illustrate sufficiently 
well that work on these materials is an interesting field of condensed matter 
physics. 
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